clear;
[y, Fs] = audioread('recoding.wav');  

Fs=48000;

n=length(y);

t=0:1/Fs:5-1/Fs;

f_noisy=3000;
noisy = 0.01*cos(2*pi*f_noisy*t');

y=y+noisy;

% 打印波形 
subplot 221;
plot(y);  
title('原始信号波形');  
xlabel('Sample');  
ylabel('Amplitude');  
ylim([-1 1]);
  
% 读取音频文件  



hw=fft(y,n);%对录制信号进行傅里叶变换
y_mag=abs(hw);
f=Fs*(0:n-1)/n;
subplot 223;
plot(f,y_mag);
xlabel('频率(hz)');title('原始信号频谱');
xlim([0 4000]);
ylim([0 250]);
  

IIR_Filter = iir();
y_f = filter(IIR_Filter,y);

subplot 222;
plot(y_f);
title('滤波后波形');  
xlabel('Sample');  
ylabel('Amplitude');  
ylim([-1 1]);

Hw_f = fft(y_f,n);
y_mag_f=abs(Hw_f);

subplot 224;
plot(f,y_mag_f);
xlabel('频率(hz)');title('滤波信号频谱');
xlim([0 4000]);
ylim([0 250]);
  
% 播放音频  
sound(y, Fs);
pause(5)
sound(y_f, Fs);

